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UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



AUXILIARY 
LIBRARY ROUTINE 



VIO - 290 



Generate a Random Normal Deviate (SADOI Only) 

26 plus 10 words of permanent storage at S3. Routine V 9 

must be located at symbolic address (V9). Routine S 5 

must be located at symbolic address (S5) 

About k5 milliseconds 



p + 1 



F5 pF 



26 qF 



q = address of this program 



Control is returned to the right hand side of p + 1 with 
X/4 in the accumulator. X is the random normal deviate, 
i.e. P[X < a] = r a -i- e" x / 2 dx. 



S 



¥2*. 



-oe 



It is suggested that the user enter routine (V9) about 
100 times before entering this routine for the first 
time , so as to discard the first 500 numbers generated 
by (V9). On the average, 6 random numbers generated by 
(V9) are used to get one random normal deviate. At the 
end of a given run using this routine the user can print 
out as sexadecimal characters (with 82 ^OF orders) the 
last numbers generated by V9. These are in S3 to 4S3. 
When this routine is next used these numbers can be input 
into S3 to 4S3 (with 8l kOF orders) and used as starting 
numbers* In this way a different sequence of numbers can 
be obtained on each run. 

Let y be a random variable which is uniformly distributed 
on the interval (-k,k). Let Z be a random variable 
which is uniformly distributed on (0,1 ). The random 
variable X = (y 
distributed for, 



e" y ' > Z) is approximately normally 



- 2 - 



2 2 / 

P[X < t J = P[y < t | e" y / 2 > Zj = p Jy< t, e" y /2 > p] 

P[e" y / 2 > Z] 



.-y 2 /2 



^1/8^ d*dy f\/8e-y% y £^' ,2 <* 



-U •'o 



-y 2 /2 



= n i 1 1 .i 



dZdy 



-4 



1/8 e-y 2 / 2 dy 






-4 






y -4 



-y 2 /2 

e y '^dy 



REFERENCE! 



Hence X Is approximately normally distributed. The 
approximation being due to the truncation of the normal 
distribution to (-4,4). m the machine a random number 
(y) is generated on (-4,4) gnd a random number (z> is 
generated on (0,1). If e " y ' 2 > z, y is accepted as a 
random normal deviatef otherwise y is rejected, a new 
y and Z are generated, and the comparison tried again. 
Actually y is accepted if y 2 + 2 £n Z < 0. The probability 

/4 p 
1/8 e" y / 2 dy =*■§£ = .3125. 

von Neumann, "Monte Carlo Methods". 
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LOCATION 


ORDER 


NOTES PAGE 1 V 10 




OOK(VIO) 









k2. 12L 
F5 3L 




Plant link 


1 


1*2 3L 
F5 6L 






2 


^2 6L 
LO 23L 






3 


32 13L 




+ -» go get 10 more random numbers 




L5 F 




n = 0, 1, 2, 3, J* 


k 


kO 24L 




xA 


5 


50 24L 

75 2^L 

. *K) 25L 




2 x 


6 


00 IF 




Waste 




L7 10S3 


by 22L, 2L 


< y < 1 


7 


50 7L 
50 7L 




Waste 


8 


26 (S5) 

00 IF 




2~ 5 £n y < 
2-^ny 


9 


36 12L 
00 IF 




+ -> accept, J2 Jcji y| > 32 

fly 


10 


36 12L 




+ -> accept |2/.n y| > 16 > x 2 




lA 25L 




2"^(x 2 + 2^n y) 

+ -» reject x > 2Jtn y - accept 


11 


32 L 






36 12L 




waste 


12 


L5 24L 




x/k in A 




22 F 




link 


13 


00 F 








50 13L 




get 5 random numbers 


Ik 


26 (V9) 
L5 S3 




random on (-1,1) 


15 


*K) 5S3 




put them in 5S3 




L5 1S3 




to 9S3 


16 


kO 6S3 








L5 2S3 






' 







LOCATION ORDER 

XT .to'7B3. 
L5 3S3 

18 kO 8S3 
L5>S3 

19 40 9S3 
50 19L 

20 26 (V9) 
L5 l^L 

21 k2 3L 
F5 18L 

22 42 6L 
22 3L 

23 00 IF 
L? 10S3 

2k Temp, store for x/k 

25 " " " x 2 /l6 



NOTES 



PAGE 2 



V 10 



accept if x + 2Ln y < 

get 5 more i.e. if x 2 < -2j£n y 
random numbers. 

xs R(-4,4) 
72 R(0,1) 



z=(x | e" x I 2 > y) 



